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Abstract.  We  consider  a  two-dimensional  convection  model  augmented  with  the  rota¬ 
tional  Coriolis  forcing,  Ut  +  U  ■  VXU  =  2kU±,  with  a  fixed  2k  being  the  inverse  Rossby 
number.  We  ask  whether  the  action  of  dispersive  rotational  forcing  alone,  C/J_,  pre¬ 
vents  the  generic  finite  time  breakdown  of  the  free  nonlinear  convection.  The  answer 
provided  in  this  work  is  a  conditional  yes.  Namely,  we  show  that  the  rotating  Euler 
equations  admit  global  smooth  solutions  for  a  subset  of  generic  initial  configurations. 

With  other  configurations,  however,  finite  time  breakdown  of  solutions  may  and  actually 
does  occur.  Thus,  global  regularity  depends  on  whether  the  initial  configuration  crosses 
an  intrinsic,  0(1)  critical  threshold,  which  is  quantified  in  terms  of  the  initial  vortic- 
ity,  w0  =  V  x  U0,  and  the  initial  spectral  gap  associated  with  the  2x2  initial  velocity 
gradient,  ??0  :=  A2(0)  —  Ai(0),  Aj(0)  =  Xj(VUo)-  Specifically,  global  regularity  of  the 
rotational  Euler  equation  is  ensured  if  and  only  if  4fcwo(a)  +  770(a)  <  4 k2,  Va  €  IR2  . 

We  also  prove  that  the  velocity  field  remains  smooth  if  and  only  if  it  is  periodic.  An 
equivalent  Lagrangian  formulation  reconfirms  the  critical  threshold  and  shows  a  global 
periodicity  of  velocity  field  as  well  as  the  associated  particle  orbits.  Moreover,  we  observe 
yet  another  remarkable  periodic  behavior  exhibited  by  the  gradient  of  the  velocity  field. 

The  spectral  dynamics  of  the  Eulerian  formulation,  [20],  reveals  that  the  vorticity  and 
the  divergence  of  the  flow  evolve  with  their  own  path-dependent  period.  We  conclude 
with  a  kinetic  formulation  of  the  rotating  Euler  equation. 
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mulation. 
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1.  Introduction  and  statement  of  main  results 

Finite-time  breakdown  is  a  familiar  trademark  of  nonlinear  convection  mechanism. 
Consider  the  canonical  example  of  an  TV-dimensional  system  of  free  transport  equations, 
Ut  +  U  ■  WXU  =  0.  It  follows  —  consult  Corollary  2.2  below,  that  the  solution  U(t ,  •)  will 
lose  its  initial  regularity  at  a  finite-time  if  and  only  if  an  eigenvalue  of  the  initial  velocity 
gradient  crosses  the  negative  real  axis,  i.e.,  if  and  only  if  there  exists  at  least  one  eigen¬ 
value,  A(0,  a;)  :=  X(9Uq^),  such  that  A(0,  a;)  G  IRC  Consequently,  finite-time  breakdown 
is  a  generic  phenomenon  for  the  free  nonlinear  transport.  Thus,  for  example,  irrotational 
initial  data  xU( 0,x)  =0  —  where  all  eigenvalues  A j(t,  x)  remain  real,  will  necessarily 
lead  to  hnite-time  breakdown,  except  for  non-generic  cases  where  Aj(0,  a;)  >  0,  Vj,  x, 
requiring,  in  particular,  that  the  initial  divergence  is  globally  positive,  •  U(Q,x)  >  0. 
This  general  TV- dimensional  scenario  is  completely  analogous  to  the  ID  inviscid  Burgers’ 
equation,  Ut  +  UUX  =  0,  where  solutions  of  the  latter  will  necessarily  reach  a  hnite-time 
breakdown  except  for  the  non-generic  case  of  monotonically  increasing  initial  data. 

Physically  relevant  models  are  governed  by  the  fundamental  Eulerian  convection  equa¬ 
tion  augmented  by  proper  forcing  F, 

(1.1)  Ut  +  U-  VXU  =  F. 

Here,  there  is  a  competition  between  the  hnite-time  breakdown  dynamics  driven  by  non¬ 
linear  convection  and  the  balancing  act  of  nonlinear  forcing,  F.  Different  models  show  up 
in  different  contexts  dictated  by  the  different  modeling  of  such  forcing.  Three  prototypes 
are  dissipation,  relaxation  and  dispersion.  It  is  well  known  that  if  (1.1)  is  augmented  with 
a  sufficiently  large  amount  of  either  dissipation  or  relaxation,  then  (1.1)  admits  a  global 
smooth  solution  for  a  rich  enough  class  of  initial  data.  In  both  cases  of  dissipation  and 
relaxation,  global  existence  is  secured  by  enforcing  a  sufficiently  large  amount  of  energy 
decay.  Dispersive  forcing,  however,  is  different.  The  dispersive  KdV  equation,  for  exam¬ 
ple,  Ut  +  UUX  =  Uxxx,  is  a  case  in  point.  It  admits  global  smooth  solution  while  keeping 
the  L2-energy  invariant  in  time. 

1.1.  The  rotational  model.  In  this  paper,  we  study  the  regularity  of  the  2D  convection 
model  augmented  by  rotational  forcing, 

(1.2)  Ut  +  U  ■  WXU  =  2UU,  .1  :=  (  ^  J  )  , 

subject  to  initial  conditions,  U(0,x)  =  U0(x).  Here  2k  =  e^1  where  e  is  the  Rossby 
number  [24],  e  =  expressed  in  terms  of  the  characteristic  length  L,  characteristic 
speed,  U,  and  the  amplitude  of  angular  velocity  D  of  the  rotating  body,  consult  [16,  23]. 
With  these  parameters  the  system  evolves  on  a  characteristic  time  scale  t  ~  L/U.  The 
system  admits  a  global  energy  invariant  in  time,  which  is  independent  of  the  amplitude 
of  rotation  encoded  by  the  constant  k  on  the  RHS  of  (1.2).  To  see  this,  we  note  that  (1.2) 
is  formally  equivalent  to  the  extended  3x3  system, 

(1.3)  dtp  -f  Vj  •  (pU)  =  0,  x  G  IR“,  t  G  IR+, 

(1.4)  dt(pU)  +  V*  •  ( pU  ®  U)  =  2kpJU, 


Rotation  prevents  finite-time  breakdown 


3 


which  are  the  usual  statements  of  conservation  of  mass  and  Newton’s  second  law,  gov¬ 
erning  the  local  density  p  =  p(t,x )  and  the  velocity  field  U  :=  (u,v)(t,x).  The  usual 
manipulation,  —  l/2\U\2x  (1.3)  +  f/Tx(1.4)  and  the  skew-symmetry  form  induced  by  the 
rotational  forcing  imply 

a<(^|c|2)  +  v,  ■  (tPu\u\2)  =  2kP(u,  ju)  =  o. 

The  global  invariance  of  the  energy  follows 

E(t)  :=  ^  J p(t,x)\U(t,x)\2dx  =  E(0). 

The  system  (1.3)-(1.4)  coincides  with  a  simplified  version  of  the  2D  shallow-water  equa¬ 
tions  (SWE),  lacking  the  additional  pressure  terms.  Although  (1.3)-(1.4)  should  not  be 
claimed  as  a  faithful  approximation  to  the  general  SWEs,  it  does  arise  as  a  meaningful 
simplification,  for  example,  when  centrifugal  forces  are  counterbalanced  by  underlying 
gravity  waves,  see  e.g.  [11,  15].  The  only  remaining  forcing  in  (1.3)-(1.4)  is  the  rotational 
Coriolis  forcing,  and  our  main  quest  in  this  paper  is  whether  the  action  of  dispersive  ro¬ 
tational  forcing  alone  prevents  the  generic  finite  time  breakdown  of  nonlinear  convection. 
The  answer  outlined  in  Section  4  is  a  conditional  yes.  Namely,  we  show  that  (1.2)  admits 
global  smooth  solutions  for  a  subset  of  generic  initial  configurations,  Uq.  With  other  ini¬ 
tial  configurations,  however,  the  finite  time  breakdown  of  solutions  may  -  and  actually 
does  occur.  Thus,  global  regularity  depends  on  whether  the  initial  configuration  crosses 
an  intrinsic,  0(1)  critical  threshold,  which  is  quantified  in  terms  of  the  initial  vorticity, 
ca0  :=  V,  x  U0  and  the  initial  spectral  gap,  T0  :=  (A2(0)  —  Ai(0))2. 

Theorem  1.1.  [Critical  threshold  for  rotation  forcing.]  Consider  the  2D  rotational  flow 
(1.2)  with  k  >  0.  Then  the  solution  of  (1.2)  with  initial  data  17(0,  x)  =  Uo(a)  remains 

x=a 

smooth  for  all  time ,  — oo  <  t  <  oo,  if  and  only  if  the  initial  data  U0  satisfy 
(1.5)  *0(q;)  :=  4 k[k  —  o;0(o;)]  —  T0(a)  >0,  V«  6  IR2. 


Moreover,  if  X(t)  =  X(t,a )  is  the  particle  path  governed  by  Xt  =  U(t,X ),  A"(0,a)  =  a, 
then  the  vorticity,  oj(t )  =  u(t,a )  and  the  divergence,  d(t )  =  d(t,a )  :=  divxU(t,a )  form  a 
periodic  orbit  in  phase  space  (u,d)  G  IR2,  with  a  period,  T  =  T( a),  given  by 

-  _ 2  r/2  de 

k  d — 7r / 2  (00  1  +  0o )  +  (0o  1  —  0o)sin0 
Here  90  =  9(a)  <  1  is  determined  by  the  initial  data 


v/1  +  8/cpo  - 1  V*o  ,  ,/ns 

9o  =  rz — : — Tyi —  ;  -  ,  Po  ■=  — - — — - ~2>  d0  =  d(0). 

\/l  +  Sfepo  +  1  dg  +  (yfiQ  —  2k) 


Several  remarks  are  in  order. 

1.  We  note  that  the  critical  threshold  (1.5)  is  independent  of  the  initial  divergence 

d0  :=  divxU0. 

2.  Let  us  point  out  that  system  (1.2)  could  be  viewed  as  a  crossroad  between  the  2D 
shallow-water  equations  and  the  so-called  2D  pressureless  equations,  e.g.  [1,  2,  3,  4,  8,  14], 
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corresponding  to  Theorem  1.1  with  k  —  0, 

(1.8)  dtp  +  V*  •  (pU)  —  0?  x  G  IR2,  t  G  IR+, 

(1.9)  dt(pU)+VX'(pU®U)  =  0. 

According  to  Theorem  1.1,  the  pressureless  system  admits  a  global  smooth  solution  for¬ 
ward  in  time,  and  respectively  —  reversible  in  time,  if  and  only  if  Aj(0)  ^  IR~,  and 
respectively  —  if  and  only  if  Aj(0)  ^  IR.  The  latter  is  equivalently  expressed  in  Theorem 
1.1  by  the  requirement  T0  :=  (A2(0)  —  Ai(0))2  <  0. 

3.  In  particular,  (1.2)  does  admit  global  smooth  solutions  with  negative  initial  diver¬ 
gence  in  contrast  to  the  free  transport  (k  =  0)  equation  discussed  in  the  introduction. 
It  follows  that  rotation  prevents  finite  time  breakdown,  either  by  a  large  Coriolis  forcing 
(k  »  1)  or  by  a  large  initial  rotation  (T0  <<  0). 


1.2.  Global  solutions  are  periodic.  We  continue  with  couple  of  remarks  on  the  peri¬ 
odicity  of  the  global  solutions  discussed  in  theorem  1.1.  If  we  set  ^-independent  initial 
data,  then  (1.2)  is  reduced  to  the  one-dimensional  system 


ut  +  uux  =  2  kv, 
vt  +  uvx  =  —2  ku: 

with  critical  threshold  (u'0(a))2  —  4kv '0(a)  <  4k2.  To  interpret  Theorem  1.1  in  this  sim- 
plihed  setting,  we  observe  that  the  gradient  (u,d)  :=  (~vx,ux)  solves  a  coupled  system 

(1.10)  (dt  +  udx)u  +  dw  =  2  kd, 

(1.11)  (dt  +  udx)d  +  d2  =  —2  ku, 


and  a  straightforward  computation  reveals  the  global  invariant  along  the  particle  path, 
X(t)  =  U(t,X),X(  0)  =  a, 


(2k  —  lu)2 
d2  +  u2 


Bq,  Bq  —  Bo(a) 


(2k  —  (jj0(a))2 

dl(a )  +  a;0(a)2’ 


The  critical  threshold  statement  in  this  case  reads  Bq  >  1,  stating  that  the  gradient 
(u,  d)  forms  a  closed  elliptical  orbit  in  the  phase  plane  (whereas  Bq  <  1  corresponds 
to  unbounded  parabolic/hyperbolic  orbits).  Following  the  analysis  in  Section  4,  we  also 
obtain  a  path-dependent  period  for  the  gradient 


t=2  r'2 _ m _  Vfii - 1 

^  ./ — ?r/2  («0  1  ■+  Do)  +  )Dq  1  —  Oo)sin6  V B o  +  1 


Such  path-dependent  period  of  the  gradient  reflects  the  fact  that  its  governing  system 
(1.10), (1.11),  is  a  nonlinear  perturbation  of  the  harmonic  oscillator.  As  the  Rossby  number 
approaches  zero,  however,  k  »  1,  90  ~  1,  and  the  above  path-dependent  period  T  is 
approaching  the  global  inertial  period  n/k  (the  harmonic  oscillator  period). 

As  we  shall  see  in  Theorem  1.2  below,  sub-critical  initial  data  yield  smooth  velocity 
fields,  U(t,X(t)),  with  time  period  T  =  n/k,  or  —  expressed  in  terms  of  the  original 
non-scaled  time  units,  a  period  T  =  T  U/L  —  7r/fh  This  period  of  the  particle  orbits  is 
related  to  the  global  revolution  of  the  plane.  Theorem  1.1  points  out  yet  another  remark¬ 
able  property  for  a  portion  of  the  gradient  of  U(t ,  •),  namely,  the  divergence  d(t,a),  and 
the  vorticity  u(t,  a)  which  exhibit  a  local  period  dictated  by  the  unique  initial  parameter, 
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8kp0.  It  is  instructive  to  compute  the  period  predicted  in  Theorem  1.1,  using  configu¬ 
rations  similar  to  those  encountered  in  various  applications.  Let  us  illustrate  a  couple 
of  examples  taken  from  [16].  For  the  Gnlf  Stream,  with  the  Rossby  number  e  =  0.07, 
L  =  100  km  and  U  =  1  m/sec,  we  find  that  the  vorticity  and  divergence  of  the  flow 
keep  repeating  themselves  every  T  =  T  L/U  ~  11.7 hrs;  for  the  weather  system  we  have 
e  =  0.14  with  L  =  1000 km,  U  =  20m/ sec,  and  the  vorticity /divergence  exhibit  a  period  of 
T  =  T  L/U  ~  12.2 hrs.  It  is  also  interesting  to  see  how  this  gradient  period  be  influenced 
by  the  small  Rossby  numbers.  After  rescaling  we  may  assume  initial  configuration  such 
that  do  ~  ouo  ~  1,  for  which  a  small  Rossby  number  yields  io  4k2,  po  r 2k  and  hence 
do  ~  1-  Restored  in  terms  of  the  original  time  scale,  T  =  T  L/U,  the  period  is  given  by 
the  do-dependent  elliptic  integral, 

T  =  2  r12 _ M _ ^  Jr 

fi  J-jr/2  1  +  Do)  +  (Do 1  —  Do)sinD 

which  is  close  to  the  inertial  period  when  the  Rossby  number  is  small.  For  the  earth  core, 
for  example,  we  have  T  =  11.95 hrs  with  e  =  2  x  10-7,  (L  =  3000 km  and  U  =  0.1  cm/ sec), 
whereas  for  Jupiter’s  Red  Spot  we  have  a  Rossby  number  e  =  0.015  (with  L  =  104km 
and  U  =  0.1  cm/ sec)  and  the  velocity  gradient  period  T  ~  5.13 hrs.  We  should  point  out 
the  difference  between  the  period  of  the  velocity  field  vs.  the  velocity  gradient  periods, 
which  is  due  to  our  tracking  of  the  flow  dynamics  along  the  particle  path. 

Of  course,  one  should  not  expect  the  current  cartoon  model  to  provide  a  faithful  descrip¬ 
tion  of  the  full  model  since  other  forces  which  are  ignored  at  this  stage,  such  as  magnetic 
forces,  pressure,  etc.,  play  a  decisive  role  in  the  dynamics  of  the  problem.  Nevertheless, 
the  ‘pure’  rotational  model  is  interesting  in  its  own  sake,  in  particular  since  the  rotational 
flow  is  predicted  to  be  periodic  once  smooth  solutions  are  secured  for  subcritical  initial 
data.  Surprisingly,  the  periods  computed  above  fall  within  the  physically  relevant  range. 
It  will  be  challenging  to  refine  the  estimates  for  these  periods  by  taking  into  account  other 
forces  which  should  complement  the  rotational  model. 

To  put  our  study  in  a  proper  perspective  we  recall  a  few  of  the  references  from  the 
considerable  amount  of  available  literature  on  the  global  behavior  of  nonlinear  convection 
driven  by  rotational  forcing.  Let  us  mention  the  rotating  Shallow- Water  (SW)  model 
studies  in  [13,  17,  25]  and  the  rotating  incompressible  Euler  and  Navier-Stokes  equations 
in  [5,  6,  10,  7].  The  common  feature  of  these  studies  is  rotation  dominated  flows  with 
sufficiently  small  Rossby  number  e.  The  flow  structure  has  been  extensively  studied  in 
terms  of  e  -C  1.  In  particular,  based  on  averaging  the  interaction  of  fast  waves  of  the 
rotating  Euler  equation,  two  dimensional  structures  were  shown  to  emerge  in  the  limit 
e  — >  0,  consult  for  example  [5,  13].  For  bounds  on  the  vertical  gradients  of  the  Lagrangian 
displacement  that  vanish  linearly  with  the  maximal  local  Rossby  number  we  refer  to  [10]. 
It  is  well  known  that  large-scale  atmospheric  (or  oceanic)  fields  are  in  permanent  process 
of  Rossby  (or  geostrophic)  adjustment,  [22],  A  nonlinear  theory  of  geostrophic  adjustment 
for  the  rotating  SW  model  for  small  Rossby  number  is  developed  in  [25].  The  analysis  of 
an  approximation  for  the  rotating  SWE  can  be  found  in  [17]. 

When  dealing  with  the  questions  of  time  regularity  for  Eulcrian  dynamics  without 
damping,  one  encounters  several  limitations  with  the  classical  stability  analysis.  Among 
other  issues,  we  mention  that 
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(i)  the  usual  stability  analysis  does  not  tell  us  how  large  the  perturbations  can  be  before 
losing  stability  -  indeed,  the  smallness  of  the  initial  perturbation  is  essential  to  make  the 
energy  method  work,  e.g.  the  3D  incompressible  Navier-Stokes  equation  [18];  in  particular 
consult  [9,  Theorem  1]  for  a  precise  statement  on  the  size  of  initial  perturbations  which 
give  rise  to  global  regularity  in  time. 

(ii)  the  steady  solution  may  be  only  conditionally  stable  due  to  the  weak  dissipation 
present  in  the  system,  say  in  the  ID  Euler- Poisson  equations  [12]. 

To  address  these  difficulties  we  advocated  in  [12,  19,  20],  a  new  notion  of  critical  thresh¬ 
old  (CT)  which  describes  conditional  stability,  where  the  answer  to  the  question  of  global 
vs  local  existence  depends  on  whether  the  initial  configuration  crosses  an  intrinsic,  0(1) 
critical  threshold.  Little  or  no  attention  has  been  paid  to  this  remarkable  phenomenon, 
and  our  goal  is  to  bridge  the  gap  of  previous  studies  on  the  behavior  of  rotational  Euler 
equations,  a  gap  between  the  regularity  of  Eulerian  solutions  in  the  small  and  their  finite¬ 
time  breakdown  in  the  large.  The  critical  threshold  (CT)  was  completely  characterized 
for  the  ID  Euler-Poisson  system  in  terms  of  the  relative  size  of  the  initial  velocity  slope 
and  the  initial  density;  consult  [19,  26]  for  the  CT  for  the  convolution  model  for  conser¬ 
vation  laws;  Moving  to  the  multi-D  setup,  one  has  first  to  identify  the  proper  quantities 
which  govern  the  critical  threshold  phenomena.  In  [21]  we  have  shown  that  these  quan¬ 
tities  depend  in  an  essential  manner  on  the  eigenvalues  of  the  velocity  gradient  matrix, 
A  (VXU). 

1.3.  On  the  Lagrangian  and  kinetic  formulations.  The  critical  threshold  for  the 
current  rotation  model  can  be  also  obtained,  in  a  straightforward  manner,  through  a 
Lagrangian  flow  formulation.  This  is  summarized  in  Theorem  1.2  below.  We  should  point 
out  that  it  was  the  spectral  dynamics  analysis  of  A  (WXU)  that  led  us  to  the  CT  formulation 
in  the  first  place,  which  in  turn  was  then  sought  within  Lagrangian  formulation.  In  Section 
5  we  prove 

Theorem  1.2.  [Flow  map  for  rotation  forcing.]  The  flow  map  associated  with  (1.2), 
Xa  :=  dXa/dt  =  U(t,Xa),  subject  to  initial  condition  Xq,(0)  =  a,  is  given  by 

Xa(t)  =  ^  J~1e2kJtUo(a)  +  a-  ^-J^Ufla). 

For  sub-critical  initial  data,  (1.5),  this  flow  map  is  invertible  and  periodic  with  an  inertial 
period  T  =  n/k.  The  velocity  field,  U(t,  Xa(t))  =  Uo(a )  +  2 kJ(Xa(t)  —  a),  shares  the 
same  inertial  period. 

At  this  point,  one  may  wonder  whether  this  inertial  period  is  none  other  than  the  planet 
rotation.  Actually  the  two  are  not  the  same;  the  rotating  planet  completes  one  revolution 
in  a  time  equal  to  yy,  while  the  period  of  the  particle  path  expressed  in  the  original  non- 
scaled  variables  is  T  =  TL/U  =  Thus,  the  particle  goes  around  its  orbit  twice  as  the 
planet  accomplishes  a  single  revolution,  which  is  consistent  with  the  observation  in  [11], 

The  periodicity  of  both  the  flow  map  and  the  “Lagrangian”  velocity  field  stated  in 
Theorem  1.2  enable  us  to  conclude  that  the  time-periodicity  of  the  “Eulerian”  velocity 
held  as  well  as  its  gradient  at  any  fixed  location  x. 
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Corollary  1.3.  For  subcritical  initial  data,  (1.5),  the  velocity  field  U(t,x)  and  its  gradient 
X7xU(t,x )  are  periodic  in  time  with  period  tt / k ,  i.e.  (U,\/xU)(t  +  n/k,x)  =  (£/,  V xU)(t,  x) 
for  any  t  >  0. 

Note  that  the  different  period  for  the  velocity  gradient  stated  in  Theorem  1.1  is  due  to 
its  expression  in  terms  of  Lagrangian  coordinate. 

Finally  we  conclude  in  Section  6  with  a  kinetic  formulation  of  the  current  rotation 
model. 

Theorem  1.4.  The  rotation  model  (1.2)  admits  for  the  following  kinetic  formulation 

dtf  +  f-  V,/  +  2kJf  •  V?/  =  -(M  —  /), 

e 

where  M^pU}(fi)  is  the  Maxwellian  given  by 

M  =  -f=e-^rr,  f  =  (6,f2)elR!, 
y-nT 

where  p  and  U  are  macroscopic  density  and  velocity  respectively ,  and  T  is  an  arbitrary 
fixed  temperature. 

In  Section  4  and  5  below,  we  quantify  the  same  critical  threshold  using  the  Eulerian  and 
Lagrangian  formulations,  and  it  would  be  of  interest  to  derive  the  same  critical  threshold 
directly  using  the  kinetic  formulation  in  Theorem  1.4. 

2.  Spectral  dynamics 

We  consider  a  general  nonlinear  transport  equation  (1.1),  Ut  +  U  ■  VXU  =  F,  and  we 
trace  the  evolution  of  VXU  in  terms  of  its  eigenvalues,  A  :=  A (X7xU)(t,x).  The  following 
result  is  in  the  heart  of  matter. 

Lemma  2.1.  [Spectral  dynamics,  [20,  Lemma  3.1]./  Let  A  :=  A (fij xU){t,x)  denote  an 
eigenvalue  of  VXU  with  corresponding  left  and  right  normalized  eigenpair,  {£,  r)  =  1. 
Then  A  is  governed  by  the  forced  Riccati  equatiori 

dtX  +  U  ■  WxX  +  A2  =  {£,  WxFr). 

As  an  immediate  corollary  we  obtain  the  precise  description  for  finite  time  breakdown 
of  free  nonlinear  transport. 

Corollary  2.2.  [Finite  time  breakdown  of  free  transport,  [20,  Lemma  4.1]./  The  free 
nonlinear  N -dimensional  transport 

(2.1)  dtU  +  U  ■  WXU  =  0,  xe\RN , 

admits  global  smooth  solution  forward  in  time,  t  >  0,  if  and  only  if  the  eigenvalues  of  its 
initial  velocity  gradient,  X  :=  X(VxU),  satisfy  A(0,  x)  £  IFC.  Likewise,  it  admits  a  globally 
smooth,  time-reversible  solution  for  —oo  <  t  <  oo  if  and  only  if  X(0,x)  ^  IR. 

For  the  proof,  we  note  that  the  eigenvalues,  governed  by  the  homogeneous  Riccati 
equations  (2.1),  propagate  along  the  particle  path  x  =  x(t,a), 

A(0,  a) 
t A (0,  a)  +  1 


A (t,  x) 
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We  note  in  passing  that  the  rotational  system  (1.2)  is  to  the  full  shallow- water  equations 
as  the  free  transport  model  (2.1)  is  to  the  full  Euler  equations.  The  existence  of  a  critical 
threshold  phenomena  associated  with  global  linear  forcing  model  was  first  identified  by  us 
[20] ,  although  the  exact  configuration  cannot  be  obtained  in  such  generality.  The  current 
paper  provides  a  precise  description  of  the  critical  threshold  for  the  2D  rotational  system 

(1.2) .  In  particular,  we  use  the  Spectral  Dynamics  Lemma  to  obtain  remarkable  explicit 
formulae  for  the  critical  threshold  surface  summarized  in  the  main  Theorems  1.1.  Taking 
the  gradient  of  the  velocity  equation  (1.2),  we  find  that  the  velocity  gradient  held,  X7xU, 
solves  the  following  matrix  equation 

(2.2)  dt(VxU)  +  U  ■  VX(VXU)  +  (V.rf/)2  =  2 kJVxU. 

Using  the  spectral  dynamics  Lemma  2.1,  we  obtain  the  spectral  dynamics  equations 

(2.3)  dtX\  +  U  ■  V*Ai  +  A2  =  2kX\(li,  Jri), 

(2.4)  dtX2  +  U  ■  XJ x X2  X2  =  2kX2(l2,  Jr2), 


where  A i,i  =  1,2  are  eigenvalues  of  the  velocity  gradient  held  X7xU  associated  with  left 
(row)  eigenvectors  and  right  (column)  eigenvectors  7 y.  Since  J  is  skew-symmetric  we 
have  Jri  =  cq/J  and  Jr2  =  a2lj .  Noting  that  l2r2  =  bn  =  1,  one  then  has  aq  =  (r2,  Jrq) 
and 

a 2  =  (ri,  Jr 2)  =  — (r 2,  Jri)  =  -an. 

Therefore,  (2.3)-(2.4)  now  read 

(2.5)  dtX\  +  U  ■  VxAi  +  Aj  =  2 A; A i  (r 2,  Jr i)(/i,  h), 

(2.6)  dtX2  +  U  ■  VxX2  +  X2  =  —2kX2(r 2,  Jr i)(h,  l2), 

from  which  we  deduce  that  the  spectral  gap  q  :=  A2  —  Ai  and  divergence  d  :=  A2  +  Ai, 
satisfy 

dt:q  +  U  ■  Vxq  +  dq  =  -2 kd(r2l  Jri){h,  l2) 

and 

,2  2 

(2.7)  dtd  +  U  ■  Vxd  +  =  -2 kq(r2,  Jn)(h,  l2). 


On  the  other  hand,  differentiation  of  (1.2)  yields  the  V a, [/-equation  (2.2),  i.e., 


(2.8)  {dt  +  U-  V x) 


(  ux 

Uy 

)  +  ( 

V  vx 

vv  , 

/  \ 

ui  +  vxu 


dV-r 


x 


dUy 

vxuy  +  vi 


=  2k 


-w. 


-U,, 


which  in  turn  -  using  the  LHS  of  (2.5)-(2.6)  to  express  A2  +  A|  =  (d2  +  q2)/2,  leads  to 


d2  4-  rr 

(2.9)  dtd  +  U  •  Vxd  +  1  =  - 2ku , 

(2.10)  dtu  +  U  ■  Vxiv  +  dto  =  2 kd. 


Equating  the  expressions  on  the  right  of  (2.9)  and  (2.7)  we  find 

-2 ku  =  -2 kq(r2l  Jri){h,l2)- 
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Thus,  the  scaled  product  of  the  eigenvectors  measures  the  ratio  of  vorticity  over  the 
spectral  gap  in  the  following  manner 


(2.11) 


U 


{r 2,  Jri)(li,l2)  =  —■ 
V 


When  the  spectral  gap  rj  shrinks  to  zero,  the  scaled  product  becomes  unbounded  due  to 
the  degeneracy  of  eigenvectors.  When  the  vorticity  c o  shrinks  to  zero,  (2.11)  recovers  the 
symmetry  of  VXU  which  is  reflected  through  the  orthogonality  of  i\  and  £2-  Equipped 
with  the  above  relations  we  come  up  with  a  closed  system  for  (u,  d,  77)  along  the  particle 
path  (here  and  below  '  =  dt  +  U  ■  Vx) 

uJ  T  du  =  2kd, 


d'  + 


d2 


T) 


r\  +  dr]  =  — 2k 


=  —2ku, 
du 


V 


Note  that  the  spectral  gap  may  become  purely  imaginary  when  eigenvalues  are  complex. 
To  avoid  the  discussion  on  the  complex  solution  of  the  above  system,  we  introduce  the 
following  real  variable 

T  :=  r f. 

Using  the  above  equations  we  have 

T'  =  2  777/  =  2d[— 2ku  —  T], 

Note  that  the  sign  of  2k  indicates  the  direction  of  the  rotational  forcing,  and  the  vorticity 
measures  the  rotation  in  the  flow.  In  order  to  combine  these  two  effects  we  introduce 
p  \=  4 k2  —  2ku,  and  thus  obtain  a  closed  system  for  W  :=  (ip,  d,  T)T 

(2.12)  p'  =  -dp, 

d2  +  T 


(2.13) 

(2.14) 


d'  = 


+  p  -  4 k2, 
4  k2-T]. 


T'  =  2  d[p 

We  shall  use  this  system  to  describe  the  dynamics  of  the  velocity  gradient  field.  Lin¬ 
earization  of  the  above  system  around  W*  =  (73*,d*,T*)T  gives  the  linear  system  W'  = 
A(W*)(W  —  W*)  with 


f  -d* 

-p* 

0 

A  = 

1 

-d* 

1 

2 

V  2d* 

2  p* 

-2d* 

The  corresponding  eigenvalues  of  A  at  critical  points  (t9*,0,T*)  are  Ai  =  0,  A2,3  = 

±y/—2p*.  The  classical  stability  analysis  based  on  linearization  is  not  sufficient  to  predict 
the  global  time  dynamics. 

3.  Material  Invariants 

It  follows  from  equations  (2.12)  and  (2.14)  that 

dp  —p 

dr  ~  2{p  —  Ak2  —  T)  ’ 


10 


H.  Liu  and  E.  Tadmor 


which  upon  integration  gives  the  first  material  invariant 


(3.1) 


2p  —  T  -  4 k- 


—  Co, 


Co  =  Co(o>) 


2p0  -I'o-  4 At 


Vo 


This  material  invariant  enables  us  to  reduce  the  full  system  (2.12)-(2.14)  to  the  following 
system 

(3.2)  p'  =  —pd, 


(3.3) 


d'  =  —~[d2  +  4k2  —  Cop2}. 


In  order  to  have  global  bounded  solution  it  is  necessary  to  assume  C0(a)  >  0,  i.e., 

(3.4)  T0  <  2p0  —  4 k2  =  4 k(k  —  u>o), 

for  otherwise,  (3.3)  will  be  majorized  by  the  Riccati  equation  d1  <  —  y  —2k2,  which  would 
lead  to  finite  time  breakdown.  As  we  shall  see  below  in  §4,  the  positivity  of  Co  (a)  is  also 
sufficient  for  global  bounded  solutions.  Another  material  invariant  is  obtained  along  the 
lines  of  [21,  Lemma  2.2]:  we  set  q  =  d2  to  find 


dq  d!  q  +  4  k2 
—  =  2d —  =  - 


dp 


Integration  yields 

/ q  +  4A;2  +  Cop2 

(3.5)  — 


<P 


p’  p 


—  TC>o,  Vo  =  V0(a)  :  — 


dl  +  4  k2  +  Copl 


\(t,Xa(t))  '  ’  po 

and  together  with  (3.1)  we  end  up  with  a  second  independent  material  invariant 

d2  —  T  +  2p\ 


{ t,xa(t )) 


=  Vo  (a) 


In  summary  we  have 

Lemma  3.1.  Let  p  :=  Ak 2  —  2 ku,  d  :=  Vx  ■  U  =  tr(WxU )  mid  V  :=  (A2  —  A1)2  be 
the  solution  of  the  dynamical  system  (2.12)-(2.1f),  associated  with  the  rotational  system 
(1.2).  Then  we  have  the  following  material  invariants  along  particle  path  ( t ,  Xa(t)), 


(3.6) 

(3.7) 


2p  —  T  —  4k2 


P 


2po(a)-T0(a)-4k2 


d2-T 


<P 


( t,xa(t )) 


(■ t,xa(t )) 


d2o( 


a 


<^o(«) 

■  r0(a) 


p0(a) 


4.  Critical  Thresholds 


As  we  observed  earlier,  the  positivity  of  condition,  C0(a)  >  0,  is  necessary  for  global 
bounded  solution,  for  otherwise 

(4.1)  d’  < -\[Ak2  +  d2}, 

which  would  imply  that  d,  and  hence  p,  become  unbounded  in  a  finite  time.  We  shall  show 
that  the  same  positivity  condition,  Co(a? )  >  0,  is  in  fact  sufficient  for  the  existence  of  global 
bounded  solution.  For  C0  >  0,  the  reduced  system  (3.2)-(3.3)  has  two  unique  equilibrium 
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points  in  the  phase  plane  (p,d)  G  IR2,  ((p±,d)  =  (±-^=,0).  The  local  behavior  of  the 
solution  depends  on  the  properties  of  these  critical  points.  We  note  that  since  p>  =  0 
is  an  invariant  set,  then  ip0(p(t)  >  0  for  all  time,  and  we  therefore  concentrate  on  the 
solution  behavior  for  tp0  >  0,  with  the  other  case  of  <^o  <  0  being  handled  similarly.  On 
the  right  plane  tp  >  0,  the  coefficient  matrix  of  linearized  system  of  (3.2)-(3.3)  around  the 
equilibrium  point  (</?+,  d)  =  (^==,0)  is 

0  -<p%  \ 

C0<P*+  0  )  ’ 


with  purely  imaginary  eigenvalues,  ±(^/CoV9+)*-  This  means  that  the  bounded  trajectory 
is  possibly  a  periodic  solution  or  limit  circle.  Observe  that  if  (tp(t),d{t))  is  a  solution  to 
(3.2),  (3.3),  so  is  ((p(—t),—d(—t)).  Such  symmetry  implies  that  (</?+,  0)  is  a  center  and 
the  trajectory  in  a  neighborhood  of  this  equilibrium  point  is  periodic. 

In  order  to  clarify  the  global  behavior  of  the  flow  around  the  center,  we  appeal  to  the 
material  invariant  (3.5)  which  we  rewrite  as 

V+M  := 

<P 

0+(-)  is  a  positive  definite  Liapunov  function  for  ip  >  0,  and  achieves  its  global  minimum, 
V+  =  0,  at  the  equilibrium  point  ( <p\  =  ^=,0).  A  family  of  closed  orbits  in  the  phase 
plane  ( p> ,  d)  can  be  expressed  as  the  level  set  curve 


V+((p,  d)  =  Const  >  0, 

since  V+  is  material  invariant  in  the  sense  that  =  0.  Similarly,  on  the  left  plane  p>  <  0, 
one  may  use  the  Liapunov  functional 


V-(<P,d)  = 


d2  +  (y/g^  +  Zkf 

-v 


whose  level  set  curve  determines  a  family  of  closed  orbit  on  the  left  plane  centered  around 
(< p>*_  =  —  A=,0).  The  global  behavior  of  the  solutions  is  summarized  in 


Lemma  4.1.  [Bounded  solutions  are  periodic.]  The  system  (2.12)-(2.14)  admits  a  global 
bounded  solution  if  and  only  if  its  initial  data  (<£>o,  do,  To)  lie  in  the  sub-critical  region  ( 
—  independent  of  do),  where 

(4.2)  T0  <  2ipo  ~  4 k2. 


Moreover,  the  bounded  solutions  of  (2. 12)- (2. 14)  are  necessarily  periodic.  The  periodic 
orbit  on  the  right  plane  ip  >  0  lies  on  the  ellipse  d2  +  (v^oV?  —  2 k)2  =  (Vq  —  4k\fC])ip, 
where  C0  and  V0  are  determined  by  the  initial  data  in  (3.1)  and  (3.5). 


Proof.  Equation  (3.3)  shows  that  if  C0  <  0  then  (4.1)  holds  and  the  divergence  blows 
up.  Thus,  the  positivity  of  C0(a )  Va  G  IR  —  namely,  the  sub-critical  condition  (4.2),  is 
necessary  for  global  regularity.  We  turn  to  show  that  it  is  sufficient.  Indeed,  if  C0  >  0 
then  the  second  material  invariant  in  (3.5)  takes  the  equivalent  form 

df  +  ( \[Cop>  —  2  k)~  =  {[Do  —  4:k\/Co)ip. 
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The  corresponding  solution,  (<p,d),  travels  along  this  elliptical  orbit  while  being  kept 
bounded.  In  fact,  ip  and  d  are  periodic  and  by  the  invariance  of  (3.6)  or  (3.7),  T  shares 
the  same  period  along  the  elliptical  orbit.  Finally,  it  remains  to  show  the  boundedness 
of  the  whole  velocity  gradient,  VXU.  The  latter  follows  from  the  boundedness  of  the 
divergence  d  along  the  lines  of  [21]  [Lemma  2.1].  To  this  end  we  note  that  the  ‘anti-trace’ 
of  (2.8),  r  :=  vx  +  uy,  satisfies  r'  +  rd  =  —2 ks,  where  the  ‘ anti- vorti city’  s  :=  ux  —  vy  is 
governed  by  s'  +  sd  =  2 kr.  Solving  the  2x2  couples  system  in  terms  of  the  divergence  d 
we  obtain 

-  =  tan(tan_1(r0/so)  —  2 kt),  r 2  +  s2  =  (r(2  +  Sg)  exp  ^—2  J  d(£)d{^J  . 

Thus,  s 2  +  r2  remain  bounded.  In  fact,  by  the  periodicity  of  d  and  its  symmetry  about 
the  axis  d  —  0  we  conclude  that  s 2  +  r2  shares  the  same  period  with  d.  A  more  precise 
statement  follows  from  the  identity  s2  +  r2  =  T  +  c o2  which  imply 

r  =  sm(tan_1(r0/so)  —  2 kt)  x  VT  +  c a2,  s  =  cos(tan_1(r0/s0)  —  2 kt)  x  VT  +  c o2. 

Being  the  product  of  two  periodic  functions  with  the  corresponding  periods  7 r/k  and  T, 
we  conclude  that  r  and  s  are  periodic  if  the  ratio  of  these  periods, 

jl=2  r _ ^ _ 

n/k  J0  (d0  +  d^1)  +  (90  -  6o1)cos(^l£), 

is  a  rational  number.  In  this  case,  the  overall  gradient  VXU  is  periodic  with  a  integer 
multiple  of  7 r/k  as  its  period.  This  completes  the  proof.  □ 

Once  we  identified  bounded  solutions  as  periodic,  the  next  step  is  to  seek  the  period  for 
each  periodic  orbit. 


Lemma  4.2.  The  period  of  each  bounded  orbit  associated  with  (2.12)-(2.1f)  is  given  by 

fAn.  —  2  r/2  dd 

k  do  +  d0  1  +  (do  1  —  do)sind 

Here  d0  =  d0(a )  <  1  is  given  by 

(4.4)  O0  :=  LSfpxl, 

where  p0  is  determined  by  the  initial  data 

/  \  a/2^o  -  r0  -  4 k2 


p0(a)  = 


d2  +  (y/2<p0-Ak2  -T0  -  2k)2 


Proof.  Due  to  the  symmetry  it  suffices  to  compute  the  half  period.  The  intersection  points 
of  the  ellipse  V (tp,  d)  —  V (<p±,  0)  with  d  =  0  can  be  written  explicitly  in  terms  of  the  initial 
data 

(4-5)  T-  =  -7^0,  T+  = 

VO  VO 

The  trajectory  from  (<p_,0)  to  ((p+,  0)  in  the  lower-half  (<p,  d)-plane  is  given  by  (3.5) 
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By  (3.2),  ip  =  —  (pd  along  this  trajectory,  whose  period  is  therefore  given  by 

rp+  dtp  2  r+  ds 


(4.6) 


T  =  2 


9_  —fd  VCo  Jip_  Sy/ {ip+  -  s)(s  -  ip-) 


_  y-+y+  I  lp+-<p- 


Let  s  =  +  2 


r;  using  the  expression  of  ip±  in  (4.5)  we  conclude 


T  —  4 


dr 


Vc,  7_i  [< p _  +  <p+  +  (< p+  -  <p-)t\V  1  -  T<2 
4  r'2  do 

VV  J-n/2  V-+V+  +  (<p+  -  <P-)sin0 


2  r/2 


do 


k  J-tt/2  @o  +  0o  1  +  (60  1  —  9o)sin6 
which  gives  the  desired  result  in  (4.3). 


□ 


5.  Flow  Map 

For  the  smooth  flow  we  may  further  study  the  structure  of  the  flow  map.  Assume 
x  =  Xa(t)  is  the  flow  map  started  at  the  initial  position  a,  then  one  has 

dX 

Xa:=^  =  U(t,Xa),  Xa(0)=a, 

and  the  momentum  equation  can  be  written  as 

Xa  =  2kJXa. 

Integration  once  gives 

Xa  =  U,{a)  +  2kJ(Xa  —  a), 

where  U,(a)  is  the  initial  velocity  at  location  a.  The  above  equation  leads  to  the  flow 
map  expression 

(5.1) 


-Mi)  =  +  a  -  ±.j-‘Uo(a). 


This  flow  map  determines  the  unique  smooth  velocity  held  if  and  only  if  the  indicator 
matrix, 

r  ■—  —  T  1  7-1/  T  2kJt\y?  TT 

w  (I~e  )VA 

remains  nonsingular.  Noting  that  J~l  =  —  J  and 

2k  it  f  cos(2kt )  sin(2kt)  \ 

6  “  '  —sin(2kt)  cos(2kt)  )  ’ 


we  have 


r  _  T  J_  j  (  1-  cos(2kt )  - sin(2kt )  .  _  „ 

+  2k  1  sin(2kt )  1  —  cos(2kt )  '  a  °' 
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Vq,£7o 


Hence,  with  U  =  ( u ,  v )  we  find 

9 hr  =  9hT  +  (  sin(2kt )  1  “  cos(2kt) 

y  cos(2kt )  —  1  sin(2kt ) 

2/c  +  v0x  +  uxsin(2kt )  —  voxcos(2kt )  uoy  +  uoysin(2kt )  —  voycos(2kt ) 

— Wox  +  v0xsin(2kt )  +  u0xcos(2kt )  2k  —  u0y  +  v0ysin(2kt )  +  u0ycos(2kt ) 

A  careful  calculation  gives  its  determinant  as 

det(2AT)  =  4A;2  —  2kuj0  +  2det(V  aU0)  +  (2/nno  —  2det(V  aU0))cos(2kt)  +  ( 2kd0)sin(2kt ). 
Thus  r(t)  remains  nonsingular  for  all  time  if  and  only  if  det(2AT)  7^  0,  i.e., 

4/c2  —  2/ano  +  2det(y  JJq) 

<£  (- aJ (2ku0  -  2det(VU0))2  +  4Fd2,  ^ (2A;cd0  -  2det(VaCo))2  +  4A;2d2")  , 


which  is  equivalent  to 

(5.2)  (4/c2  —  2/ano  +  2det(V aU0))2  >  (2kcu0  ~  2det(V  aU0))2  +  Ak2d2Q. 

We  now  invoke  the  relation  between  the  spectral  gap  r0,  and  the  corresponding  the 
determinant  and  divergence, 

T0  =  d2  -  Adetiy  JJq). 

Applied  to  the  above  inequality  (5.2),  this  yields 

r0  <  Ak 2  —  4/ano, 


which  is  exactly  the  critical  threshold  (1.5)  stated  in  Theorem  1.1. 


Shifting  our  attention  to  the  Eulerian  framework,  we  conclude  this  section  with  a  discus¬ 
sion  on  the  time-periodicity  of  sub-critical  velocity  field  for  fixed  location  x.  The  solution 
determined  by 

U(t,  x )  =  e2kJtU0(a),  x  =  a  +  ^  J~\e2kjt  -  1  )U0(a), 

is  implicit.  For  sub-critical  initial  data  one  could  find  a  in  terms  of  (t,  x)  and  substitute 
back  to  obtain  the  velocity.  In  fact  a  is  given  by 

a  =  x  +  7 —J~1(e~2kjt  —  1  )U(t,  x ), 

ZrC 

which  gives  the  velocity  field  implicitly  determined  by 

U(t,x)  =  e2kJtu0  (x  +  ^(e~2kJt  - 1  MM))  • 

Combined  with  the  critical  threshold  condition  (1.5),  this  shows  that  for  any  fixed  x , 
U (t,  x)  and  hence  VXU ( t ,  x)  are  periodic  in  time  with  period  n/k  and  corollary  1.3  follows. 


Rotation  prevents  finite-time  breakdown 
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6.  Kinetic  formulation 

This  section  describes  a  kinetic  formulation  for  the  rotational  model  (1.2)  in  terms  of 
a  density  function,  /  =  f(t,x,£)  governed  by  the  BGK  model 


(6.1) 


a,  /  +  E-V./  +  2fcfL  •  V{/  =  -(M  -  /),  y  := 


Here  M  =  Afnm(^)  is  the  Maxwellian  given  by 


M  =  -(Le-I«-U|2/T 
\fnT 


€  =  (6,6)  e  IR2. 


The  fixed  temperature  T,  plays  no  role  in  this  pressureless  model.  The  connection  between 
the  distribution  function  /  and  macroscopic  flow  variable  is  realized  in  terms  of  the  usual 
moments  of  density  p,  momentum  m  =  pU  and  total  energy  K1  =  p|t/|2/2, 

(p.pU.E)1  =  J  VT(OM,  VR)  :=  (u.^)! 

The  conservation  principle  for  mass,  momentum  and  energy  during  the  course  of  particle 
collisions  requires  the  equilibrium  to  satisfy  the  compatibility  condition 

J (M  -  f)ij>T (£)d£  =  0, 

while  the  rotational  forcing  is  introduced  through  the  potential 

J  2k£±  ■  Vtfi/>T{t)d£  =  (0,  —2kpJU,  0). 

Indeed,  a  straightforward  computation  yields  f  2 k^T  ■  =  0;  for  the  momentum 

equation  we  compute 

[  2ke  ■  Viftdt  —  2k  [  e(6%/  -  ZidtofW 


=  2k 
=  2k 


66%/ 

-66%/ 


d{ 


—pv 
pu 

=  —2  kpJU, 

while  the  presence  of  this  forcing  does  not  change  the  energy  equation  since 

2  kje-  Vzm2dt  =  o. 

The  first  three  moments  of  (6.1)  then  yield  the  equivalent  extended  system  of  (1.3), (1.4), 


d. 


+  Va 


0 

=  (  2kpJU 
0 


The  corresponding  macroscopic  fluxes  are 


(Fp,Fm,FE)T  :=  I  iP'imfdt, 
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and  under  the  closure  /  =  M^p  Uj  we  conclude 

T  /  |t/|2X  T 

(Fp,  Fm,  FEy  =  (pU,  pU  ®  U,  y-j-)  . 
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